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ABSTRACT 



Context. Do extrasolar planets affect the activity of their host stars? Indications for chromospheric activity enhancement have been 
found for a handful of targets, but in the X-ray regime, conclusive observational evidence is still missing. 

Aims. We want to establish a sound observational basis to confirm or reject major effects of Star-Planet Interactions (SPI) in stellar 
X-ray emissions. 

Methods. We therefore conduct a statistical analysis of stellar X-ray activity of all known planet-bearing stars within 30 pc distance 
for dependencies on planetary parameters such as mass and semimajor axis. 

Results. In our sample, there are no significant correlations of X-ray luminosity or the activity indicator LxlLtoi with planetary 
parameters which cannot be explained by selection effects. 

Conclusions. Coronal SPI seems to be a phenomenon which might only manifest itself as a strong effect for a few individual targets, 
but not to have a major effect on planet-bearing stars in general. 
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1. Introduction 

The detection of extrasolar planets is one of the outstanding 
achievements in astronomy during the last 20 years. The first de- 
tected exoplanet rev ealed propert ies which were very surprising 
at that time: 51 Peg dMavor & OueloZi.1995) hosts a Jupiter-like 
planet at a distance of only 0.05 AU, thus the planet orbits its 
host star in less than five days. Since then, more than 400 other 
exoplanets have been found at the time of writing (see for exam- 
ple the Exoplanet Database at www.exoplanet.eu), both in very 
close orbits and in such more familiar from our own Solar sys- 
tem. 

With the existence of extrasolar planets established, the 
question arises what the environmental properties of such plan- 
ets may be and if they might even allow the existence of life. 
The physical properties of planets, especially in close orbits, are 
crucially determined by the irradiation from their host stars. The 
evaporation rate of a planetary atmosphere depends on its ex- 
ospheric temp erature Too , i.e., the regions where particles can 
escape freely jLammer et alj|2003h . Thus, the host star's EUV 
and X-ray radiation is the key property determining a plan- 
ets exospheric temperature. Evaporation of the planetary atmo- 
sphere has been observed for the transiting planet HD 209458b 
dVidal-Madiar et al]l2003h : the planet loses hydrogen which is 
observable in absorption spectra during the transit. 

At very close distances, one might expect also planets to in- 
fluence their host stars, in analogy to binary stars which show a 
higher activity level compared to single stars. Two different pro- 
cesses for Star-Planet-Interaction (SPI) have been put forward 
(|Cuntz et al. 2000). Planets can induce tidal bulges on the star 
with an interaction strength depending on the planetary semima- 
jor axis (oc cr^j ), which might lead to enhanced coronal activity 
via increased turbulence in the photosphere. Planetary magnetic 
fields can also interact with the stellar magnetic field (oc a'^) and 
might also induce enhanced activity via Jupiter-Io-like interac- 



tion, i.e. flux tubes which connect star and planet and heat up 
their footpoints on the stellar surface, or magnetic reconnection. 
Some observational campaigns have been conducted to investi- 
gate the existence of possible SPI: Shkolniket al. (2005) moni- 
tored the chromospheric activity of 1 3 stars via Ca n H and K line 
fluxes and found indications for cyclic activity enhancements in 
phase with the planetary orbit for two of these stars. The activ- 
ity enhancements appeared once per planetary orbit, suggesting 
magnetic instead of tidal interaction. Howeve r, measurements 
obtained three years later dShkolnik et al.ll2008h showed that the 
activity enhancements had switched to a cycle in phase with the 
stellar rotation period instead. 

The coronal activity of planet-b earing stars has been inves- 
tigated in a first systematic study bv lKashvap et al.l (2008). The 
authors claim an over-activity of planet-bearing stars of a factor 
of four compared to stars without planets, but their study had to 
include upper limits for a large number of stars since less than 
one third of the stars in their original sample were detected in 
X-rays at that time. A dedicated campaign to search for mag- 
netic SPI in the case of HD 179949, one of the stars which 
IShkolnik et alj (l2005h found to have cyclic activ i ty cha nges in 
the chromosphere, was conducted by ISaar et aO (l2008h . These 
authors found spectral and temporal variability phased with the 
planetary orbit, but some of that might also be induced by intrin- 
sic stellar activity variations, since the stellar rotation period is 
poorly known (P, = 7 - 10 d). 

Up to now, the observational basis of stellar coronal activity 
enhancements due to close-in planets is not sound enough to es- 
tablish or reject the possibility of coronal SPI. In order to adress 
this issue we conducted an X-ray study of all planet-hosting 
stars within a distance of 30 pc with XMM-Newton which have 
not been studied with ROSAT before. In this fashion a volume- 
limited complete stellar sample can be constructed. 
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2. Observations and data analysis 

As of December 1st 2009, a total of 72 stars within 30 pc 
distance have been detected which are known to to harbor 
one or more planets. For some of these. X-ray properties are 
known from previous ROSAT or XMM-Newton observations, 
but for a large number of these stars X-ray characteristics were 
not or only poorly known. Therefore we observed a total of 
20 planet-hosting stars with XMM-Newton between May 2008 
and April 2009 to determine X-ray luminosities for stars which 
had not been detected before in other X-ray missions, and to de- 
rive coronal properties from spectra recorded with EPIC (both 
MOS and PN CCD detectors) especially for stars with close-in 
planets. We reduced the data with SAS version 8.0, using stan- 
dard criteria for filtering the data. We extracted counts from the 
expected source regions with radii between 10" and 30", depend- 
ing on the source signal, background conditions and the presence 
of other nearby sources. Background counts were extracted from 
much larger, source-free areas on the same chip for the MOS de- 
tectors and at comparable distances from the horizontal chip axis 
for the PN detector 

For hitherto undetected stars, which showed only a weak 
source signal in our observations, we used the source detection 
package "edetect-chain" of SAS v8.0. As stars with low X-ray 
luminosities in general have lower coronal temperatures and thus 
softer spectra, we used energy bands of 0.2-1 keV in PN and 
0.15-1 keV in MOS and merged all EPIC detectors for source 
detection. 

For the subsequent analysis of all stars we use the four en- 
ergy bands 0.2-0.45 keV, 0.45-0.75 keV, 0.75-2 keV and 2-5 ke V, 
because not all of our sample stars were detected with sufficient 
signal to noise ratio to allow spectral fitting. With the four en- 
ergy bands, we can calculate the stellar fluxes via ECFs (energy 
conversion factors) for each band more accurately than by just 
assuming a single ECF for all counts. Above 5 keV, there is very 
little to no signal present in comparison to softer energies for all 
of our stars. We calculated these ECFs by simulating spectra in 
Xspec vl2.5 for diff'erent coronal temperatures with the respec- 
tive instrumental responses and efifective areas of the detectors 
folded in. This yields reliable ECFs which vary about 25% for 
coronal temperatures above 1 MK for thin and medium filters. 
For the thick filter, the small effective area below 350 eV intro- 
duces larger errors in the ECFs already for temperatures below 
logT[K] = 6.2. 

For the error estimate on our derived luminosities, we use 
Poissonian errors on the total number of source counts, and an 
additional error of 30% to account for uncertainties in the ECFs 
and stellar variability. For stars which were not observed with 
XMM-N ewton , we u se the published X-ray luminosities from 
[Kashva p et al.l (l2008h and a dd an extra error of 4 0% on top of 
their Poissonian errors, since lKashvap et all (l2008h used a single 
ECF for their flux calculations. 

When comparing X-ray luminosities derived from XMM- 
Newton and ROSAT observations, one has to take into ac- 
count the dififerent energy bands accessible to the detectors 
(0.2-12 keV for XMM-Newton, 0.1-2.4 keV for ROSAT). For 
coronal temperatures between logT = 6.2 and 7.0, the flux 
in the ROSAT band is larger by a factor of 1.1 compared to 
the XMM-Newton band. For lower coronal temperatures be- 
tween logT = 6.0-6.2 and log T = 5.8-6.0, the flux in the 
ROSAT band is larger by a factor of a 1.5 and * 4.0, respec- 
tively, rising steeply towards even lower temperatures, since the 
spectrum shifts to energies which are inaccessible to XMM- 
Newton. With these factors, we can transform the XMM-Newton 



fluxes to the ROSAT band, identifying stars with coronal tem- 
peratures below logT = 6.2 and 6.0 by a hardness ratio of 
HR = (H -S)/(H + S)<+0A9 and -0.34 respectively, with H 
and S being the source counts in the energy bands 450 - 750 eV 
and 200 - 450 eV, respectively. 

We did not exclude flaring periods of individual stars when 
doing comparisons of X-ray luminosities or activity indicators, 
since we cannot identify flares in stars which are barely de- 
tectable and do not allow lightcurve analysis. We do, however, 
distinguish between flaring and quiescent phases for spectral 
analyses of individual stars. 

When conducting intra-sample comparisons, we will use 
only detections for Kolmogorov-Smirnov and correlation tests, 
but we will include upper limits when doing linear regressions 
of X-ray luminosities or activity indicators over planetary prop- 
erties. 



3. Sample properties 

Now we characterize our sample of planet-bearing stars within 
30 pc with respect to X-ray detection rates and X-ray surface 
fluxes of ROSAT- and XMM-Newton-detected stars as well as 
in comparison with field stars. Tables l4l and 16. 1 1 list stellar and 
planetary parameters as well as X-ray properties of the sample 
stars which have been detected with XMM-Newton and ROSAT, 
respectively. 

3. 1 . X-ray detection rate 

In total, 72 stars planet-bearing stars have been detected within a 
distance of 30 pc. 36 of these were observed with XMM-Newton 
over the last years, yielding 32 X-ray detections. For 24 addi- 
tional stars, which have not been observed with XMM-Newton, 
X-ray luminosities are known from ROSAT observations. This 
yields 56 stars with known Lx out of the total sample of 72 stars. 
In our further sample analysis, we will leave three detected stars 
out of our analysis, namely y Cep, Fomalhaut and /3 Pic; the for- 
mer being a spectroscopic binary which cannot be resolved in 
X-rays, the two latter being A-type stars for which the produc- 
tion process of X-ray emission is supposedly very different from 
later-type stars with a corona and therefore also any planetary in- 
fluence on X-ray properties should be determined by a different 
mechanism compared to stars of spectral type F and later. 

The stars within 30 pc around which planets have been de- 
tected are mainly of spectral type G or later. Fig.[T]gives the rate 
of X-ray detections versus spectral type, being 75% for F stars, 
> 65% for G stars and > 85% for K and M stars. 

3.2. X-ray surface flux 

For a subsample of our stars excluding giants, we examined 
the X-ray surface flux. The lowest flux level of XMM-Newton- 
detected stars seems to be systematically lower than for ROSAT- 
detected stars. This is not surprising, since both X-ray telescopes 
have diff'erent accessible energy bands (0.2-12.0 keV for XMM- 
Newton, 0.1-2.4 keV for ROSAT) and the integrated X-ray flux 
depends strongly on the lower energy cutoff especially for cool 
coronae as in our stars. 

Gauged to the same energy band (as described in section|2]l, 
both XMM-Newton- and ROSAT-detected stars show a limit- 
ing X-ray surface flux level near log Fsurj [ergs"' cm"^] > 4.0 
(see Fig. |2]i. For the calculation of Fj,„/, we use the stel- 
lar radii given in the exoplanet.eu database. If we compare 
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Table 1. Correlation of X-ray luminosity and Lx/Lb„i with plan- 
etary parameters; X/R; XMM-Newton/ROSAT detections. 



Fig. 1. Spectral types of planet-bearing stars within 30 pc (solid); 
X-ray detections marked as dotted lines. 
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Fig. 3. X-ray luminosity of F and G type stars. Planet-bearing 

1 .0 stars: squares, stars without planets: crosses. 



Fig. 2. X-ray surface flux of planet-bearing stars vs. planetary 
distance. XMM-Newton fluxes are shown as triangles, ROSAT 
fluxes as squares. XMM-Newton fluxes scaled to the ROSAT 
energy band; the flux level of a solar coronal hole (log Fsurf ~ 4) 
is indicated by the dotted line. 

the XMM-Newton and ROSAT surface flux sample with the 
Kolmogorov-Smirnov test, we find that both populations are 
significantly different. This is due to the selection effect that 
we proposed planet-bearing stars which were previously unde- 
tected with ROSAT (and therefore have low X-ray luminosi- 
ties) for detection pointings with XMM-Newton. This leads to 
a higher concentration of stars near the limiting surface flux 
level of log/^TO,./ [erg s"' cm"^] a 4. For the XMM-Newton 
and ROSAT subsamples of stars with a surface flux above 
logFjur/ [erg s"' cm"^] > 4.5, we find that these populations 
are statistically indistinguishable (probability for both samples 
stemming from the same distribution 71%). 

3.3. Comparison with field stars 

To check for systematic differences, we compare our sample of 
planet-bearing stars with a sample of field stars of spectral type 
F and G as available from Schmitt ( 1997) from ROSAT obser- 
vations. In Fig. [3] we show the X-ray luminosities of these stars 



over B - y of both the planet-bearing and non-planet-bearing 
sample. A Kolmogorov-Smirnov test yields that the probability 
that both samples are drawn from the same parent distribution is 
74%. The values of the activity indicator Lx/Lboi yield a prob- 
ability of 23% to be from the same distribution; this can be ex- 
plained by the fact that stars of low activity are generally chosen 
for planet search programs. 

4. Star-Planet Interactions 

Now we investigate our sample in detail for possible correlations 
of X-ray properties with planetary parameters. 

One expects that possible effects which giant planets might 
have on their host stars will strongly increase with decreasing or- 
bital distance. Also, tidal as well as magnetic intercations should 
increase with the exoplanet's mass, assuming that larger exo- 
planets are capable of producing a stronger planetary magnetic 
field. Note that closer-in planets may rotate more slowly since 
they synchronize with their orbi t, weak ening their ability to gen- 
erate magnetic fields dGrieBmeier et al.ll2004l) : however, the de- 
tails of planetary dynamos are not fully understood. 

The most interesting quantity with regards to SPI in the X- 
ray regime is the activity indicator Lx/Li,,,!. The X-ray luminos- 
ity alone varies with stellar radius independently of the activity 
level, but Lx/Lhoi is independent of such radius-induced effects. 
Any planet-induced activity changes should therefore be evident 



4 



K. Poppenhaeger et al.: Coronal properties of planet-bearing stars 



in LxlLhof, a planet-induced variation in Lx which would leave 
the ratio unchanged is rather unphysical, since Lx/Lhoi has typi- 
cal values of 10"^ for our stars. A change in Li,oi would therefore 
need 10* times more energy than the X-ray variation alone. 

We study both the X-ray luminosity Lx as well as the activ- 
ity indicator Lx/Lhoi for correlations with the innermost planet's 
semimajor axis and mass. In Table[T]we give Spearman's p rank 
correlation coefficient for various combinations. A value of 1 
(-1) means a perfect correlation (anticorrelation), means no 
correlation. The corresponding p-value gives the probability that 
the observed value of p can be obtained by statistical fluctua- 
tions. 

For the correlation analysis of Lx/Lhoi, we exclude giants 
from our sample (HD 27442 and HD 62509), since they have 
very low Lx/Lhoi values due to their optical brightness. As well 
as for XMM-Newton detections alone as for XMM-Newton and 
ROSAT detection combined, we find no correlation of the semi- 
major axis with the stellar X-ray luminosity. 

We find two possible correlations here: one of planetary mass 
with Lx and a stronger one for ajj x Mpi with Lx- Stars with gi- 
ant and close-in planets have higher X-ray luminosities than stars 
with small far-out planets. For Lx/Lt,,,!- there is a correlation with 
a~l X Mpi present in the sample of XMM-Newton detections, but 
not in the larger sample of ROSAT and XMM-Newton detec- 
tions, pointing towards the possibility that this correlation might 
be a statistical fluctuation. The probable reason for strong corre- 
lations in Lx, but weaker or absent ones in LxlLhoi is that there 
is also a strong (> 2cr) correlation between Mpi and Lj,,,!'. Stars 
with larger L/,,,; are more massive, and around massive stars, gi- 
ant planets are detected much more easily compared to small 
ones. Both correlations of planetary mass with Lx and also Lhoi 
seem to cancel out in LxlL^oi. 

Another significant correlation worth mentioning exists be- 
tween the planetary mass and the spectral type of the host star: 
small planets are prone to be found around stars of later types. 
This is basically the same trend we see between Mpi and Lboi, 
since small planets are more easily detected around low-mass 
and therefore late-type stars. 

To visualize these (non-)correlations, we perform linear 
regressions of logL^ and log(Lx/Lhoi) with either loga^/ or 
log(fl^^' X Mpi) by using the "linmix_err" routine implemented 
in IDL. As already seen in the correlation analysis, we obtain 
slopes which are compatible with zero at 1 cr level for all of the 
pictured cases except for log Lx over log(fl^^' x Mpi) (see Fig.|5l). 

Independently of any linear trend, we can test if the Lx val- 
ues of stars with close-in and far-out planets stem from the same 
distribution. Fig. |4] shows the logarithmic X-ray luminosities of 
stars with planets within Upi < 0.2 AU and stars with planets 
beyond Opi > 0.5 AU. The means of both distributions are very 
similar and not distinguishable at Icr level: < log Lx dose-in >- 
27.52 + 0.72 erg s^' and < log Lx far-out >= 27.70 + 0.80 erg s^'. 
A Kolmogorov-Smirnov test yields a probability of 84% for both 
samples being from the same distribution. However, the com- 
parison of stars with close-in, heavy planets compared to far-out 
light ones yields that the probability for both samples to have the 
same parent distribution is very small with < 1%; the average 
X-ray luminosity is higher for stars with close-in, heavy planets 
with < log Lx close-in heavy >- 27.91 ± 0.76 erg s"' compared to 
< log Lx far-out light >- 27.41 ± 0.73 ergs"', but the means are 
compatible within Icr errors. 
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Fig. 4. Histograms of X-ray luminosities for X-ray detected stars 
within 30 pc with close-in (api < 0.2 AU, upper panel) and far- 
out (api > 0.5 AU, lower panel) planets. Mean log Lx values are 
indicated by dotted lines for both samples. 

5. Properties of individual targets 

In the following, we give a short overview on our newly observed 
stars and their spectral properties. The stellar, planetary and X- 
ray properties of all planet-bearing stars which were observed 
with XMM-Newton are listed in Table g] 

5. 1 . Individual targets 

Between May 2008 and April 2009 we observed a total of 20 
planet-bearing stars in X-rays. One of these stars, SCR 1845, 
turned out to harbor a brown dwarf and not a planet; the X-ray 
characteristic s of SCR 18 45 are discussed in a separate publica- 
tion ( Robrade et al.l2010 ). The X-ray properties of the remaining 
19 stars are described briefly now. 

GJ 674, GL 86, GL 876, HD 102195, GJ 317, 55 Cnc and 
HD 99492 yielded sufficient signal-to-noise ratio for spectral fit- 
ting of the obtained EPIC spectra. They are all characterized by 
coronae with cool to moderate temperatures (details listed in sec- 
tion |52]l. GJ 674 shows one large and several smaller flares on 
timescales of ca. 5 ks. Also GL 876 shows several short flares. 
The other stars show some variability around 15-25% level. As 
an example, we show the EPIC spectrum and the corresponding 
two-temperature fit of GL 86 in Figure |6] 
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Fig. 5. X-ray luminosity and activity indicator \og{Lx I Lhoi) as a function of log Upi and log(a^/ x Mpi), respectively. XMM-Newton 
detections: triangles, ROSAT detections: squares. 
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Fig. 6. Typical EPIC spectra (upper: PN, lower: MOS) of a 
planet-bearing star; here shown: GL 86. 



HD 154345, HD 160691, HD 4308, HD 52265, HD 93083, 
51 Peg, HD 27442, HD 114386 and HD 114783 were detected in 
our exposures, but did not yield enough photons for spectral fit- 
ting. Where meaningful hardness ratios could be calculated from 



the numbers of counts, the stars proved to be soft X-ray sources, 
as one expects for nearby stars with low X-ray luminosities. Two 
of these targets show interesting characteristics: 

HD 99492 is part of a binary system, consisting of a K2 
dwarf, which is the planet-bearing star HD 99492, and a KO sub- 
giant, HD 99491. Previously, HD 99492 was assigned an X-ray 
luminosity of 27.56 erg s"' (Kashvap et al. 2008), since the dou- 
ble system was unresolved in the corresponding ROSAT point- 
ing. Our XMM-Newton pointing shows that HD 99492 is actu- 
ally the X-ray fainter part of the pair with an X-ray luminosity 
of only 26.93 ergs '. 

51 Peg is, despite a moderate X-ray luminosity of 
26.8 ergs ' which it exhibited in a ROSAT pointing from 
1992, barely detected in a deep XMM-Newton observation. 
The XMM-Newton source photons are extremely soft, which 
explains its better visibility in ROSAT and additional recent 
Chandra data, since these instruments have a larger effecti ve area 
at very low X-ray energies. Detailed analysis (Poppenhage r et all 
|2009) showed that the star is possibly in a Maunder minimum 
state. 

The stars 16 Cyg B, HD 111232, HD 217107 and HD 164922 
could not be detected in X-rays in our exposures. The upper lim- 
its for these stars wer e calculated fo r a confidence level of 99%, 
following the lines of lAvresI (|2004|) and point also towards low 
activity levels for these targets. 
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Parameter GJ 674 quiescence GJ 674 flare 



Ti (keV) 


0.14 


0.32 




U. /i 


i.j / 


Ti (keV) 


0.40 


0.81 


EM2 


2.52 


1.61 





0.41 + 0.07 


0.38 + 0.05 


Ne 


0.66 + 0.13 


0.38 ±0.18 


Mg 


0.34 ±0.14 


0.51 ±0.25 


Fe 


0.29 ± 0.05 


0.37 ± 0.12 


Xl, (d.o.f.) 


1.12 (204) 


1.00(197) 


log Lx (0.2-2.0 keV) 


27.5 


27.7 



Table 3. Spectral modeling results for GJ 674 derived from EPIC 
data; emission measure given in units of 10^" cm"^. 



5.2. Spectral properties 

Seven of the stars we observed yielded sufficient signal to noise 
ratio for spectral fitting of tfieir XMM-Newton EPIC data. Tfie 
spectra of six stars can be adequately described by a thermal 
plasma with two temperature components and solar abundances; 
the spectral fitting was performed with Xspec vl2.5 and apec 
models. The derived spectral properties are listed in Table |2] 
They are mostly dominated by very cool plasma (T ^ I MK) 
with small contributions from hotter plasma, only HD 102195 
and GL 86 have stronger contributions from a hotter component 
around 4-5 MK. 

The spectrum of GJ 674 cannot be fitted satisfactorily when 
assuming solar abundances, therefore we fit its EPIC spectra 
with two-temperature vapec models. GJ 674 exhibits a flare dur- 
ing our observation, so we conducted the spectral analysis for 
the time interval of the flare as well as for quiescent times. The 
results are given in Table[3] The temperature of both components 
rises considerably during the flare, as well as the total emission 
measure. Coronal abundances of the given elements are below 
solar values, consistent with the low photospherical iron abun- 
dance [FelH] = -0.28. 

The X-ray luminosities derived from spectral modelling fit 
the ECF-derived results in Table |4] well, justifying our error es- 
timate of » 30%, which we assumed in addition to statistical 
errors. 



6. Discussion 

6.1. Interaction or selection? 

In our data, we do not see any significant trend of the activity 
indicator LxlL^oi with the planetary semimajor axis, mass or a 
combination of both, in contrast to recent studies (iKashvap et al.l 
I2OO8 ). The only significant trend, as shown in section|4l is seen 
in the X-ray luminosity which is higher for stars with heavy 
close-in planets. Trying to explain this trend in Lx without an 
accompanying trend in LxlLtoi by SPI is problematic. The sam- 
ple stars have typical LxjLhoi values of 10"^. If planets induced 
higher Lx levels, but constant LxjLhoi ratios, the amount of en- 
ergy introduced by SPI would have to be 10^ times higher in 
Lboi than in Lx. The Lx excess of stars with close-in heavy plan- 
ets compared to stars with far-out light ones is of the order of 
5 X 10^^ erg s"' . Thus, the energy excess in Ltoi would have to be 
~ lO^'-' ergs"'. Comparing this to the typical orbital energy of a 
Hot Jupiter (~ 10"*"* erg), this would lead to obviously unpyhsical 
timescales for the planet's orbital decay of only several thousand 
years. 
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Fig. 7. Linear regression of log Lx over log(fl^^' x Mpi) for close- 
in heavy planets (solid) and far-out light planets (dotted). Both 
regressions overlap at Icr level. 

However, there is also the possibility that the trend in Lx 
is caused by selection effects: All but three of the planets in our 
sample have been detected with the radial velocity (RV) method. 
Stellar activity can mask the RV signal. Since the RV signal is 
strongest for close-in, heavy planets, we have a selection effect 
which favors detection of such planets around active stars. The 
key question is, do the data show an additional trend of Lx with 
a^' and Mpi which is not induced by the selection effect and 
could be attributed to SPI? 

To investigate this, we conduct two separate regression anal- 
yses on \ogLx over Iog(a^^' x Mpi) for close-in heavy planets 
and far-out light planets, respectively. We define close-in heavy 
planets by log(fl^^' x Mpi) > 0.5 (corresponding to a Jupiter-like 
planet at a maximum orbital distance of ^ 0.3 AU, for example) 
and far-out light planets by log(a^' x Mpi) < (corresponding to 
a Jupiter-like planet at 1 AU or a Saturn-like planet at 0.3 AU). 
As shown in Fig.|7] both regressions overlap well at Icr level, 
indicating that there is no additional activity enhancement effect 
measurable in this sample for close-in heavy planets other than 
the selection trend which also manifests itself in the subsample 
with far-out light planets. 

6.2. Is there evidence for coronal SPI? 

There are two different sc enarios for SPI: ti dal and magnetic in- 
teraction (see for exampIe lCuntz et alJ (12000) ). Tidal interaction 
will affect motions in the stellar convection zone as well as the 
flow of plasma in the outer atmospheric layers. If stellar rota- 
tion and the planetary orbital motion are not synchronous, tidal 
bulges should rise and subside on the star, causing additional 
turbulence at the footpoints of magnetic loops, leading to higher 
flaring rates, or causing outer layers of the star to corotate with 
the planet, which might enhance the stellar dynamo if Port > P*- 
Magnetic interaction is thought to be able to enhance the stellar 
activity in several ways: if planets are close enough to their host 
stars to be located inside the star's Alfven radius, a Jupiter-Io- 
like interaction can form where the planet is connected with the 
star via fluxtubes which heat the stellar atmosphere at their foot- 
points. Alternatively, magnetic reconnection events of the stellar 
and planetary magnetic field lines might supply additional en- 
ergy to the stellar atmosphere. Also the mere presence of the 
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Parameter 


GL 86 


GL876 


HD 102195 


GJ317 


55 Cnc 


HD 99492 


Ti (keV) 


0.11 


0.07 


0.09 


0.11 


0.09 


0.09 




1.63 


0.42 


8.54 


0.19 


1.39 


1.02 


n (keV) 


0.33 


0.39 


0.45 


0.50 


0.48 


0.39 




1.27 


0.06 


8.07 


0.06 


0.19 


0.20 


X]^j (d.o.f.) 


1.02(120) 


1.44(85) 


0.95 (152) 


0.77 (22) 


1.2(11) 


0.99 (8) 


log Lx (0.2-2.0 keV) 


27.6 


26.4 


28.3 


26.5 


27.1 


27.0 



Table 2. Spectral modeling results derived from EPIC data; emission measure given in units of 10^ cm 



Star 


Type 


Dist. 


my 


B - V 


Pt 


[Fe/H] 


"pi 




GTI 


net countrate" 


Fx 


logLx 








(pc) 






(d) 




(AU) 


(Mj) 


(s) 


(cts/ks) 


(erg s"' cm"-) 


0.2-2 keV 


e Eri 


K2.0 V 


3.2 


3.73 


0.88 


11.2 


-0.10 


3.39 


1.55 


10385 


5611.5± 23.3" 


1.38E-11 


28.22± 0.12 


-4.88 


GJ674 


M2.5 


4.5 


9.38 


1.53 


35.0 


-0.28 


0.04 


0.04 


15183 


1102.6± 8.6 


2.16E-12 


27.73±0.12 


-3.80 


GL876 


M4.0 V 


4.7 


10.17 


1.67 


41.0 


-0.12 


0.02 


0.02 


23436 


48.9± 1.5 


1.13E-13 


26.48± 0.13 


-5.11 


VB 10 


M8.0 V 


5.8 


9.91 


0.00 


- 


0.00 


0.36 


6.40 


10810 


23.2± 1.5 


4.91E-14 


26.30± 0.14 


-1.79 


GJ317 


M3.5 


9.2 


12.00 


1.53 


- 


-0.23 


0.95 


1.20 


11245 


14.3± 1.5 


3.30E-14 


26.52± 0.14 


-4.57 


HD 62509 


G5.0 IV-V 


10.3 


1.15 


1.00 


130.0 


0.19 


1.69 


2.90 


28759 


31.2± 1.1 


l.OlE-13 


27.11± 0.13 


-6.63 


GL86 


Kl.O V 


11.0 


7.40 


0.77 


31.0 


-0.24 


0.11 


4.01 


13071 


116.8± 3.1 


2.94E-13 


27.63± 0.12 


-5.04 


55 Cnc 


G8.0 V 


13.0 


5.95 


0.87 


42.7 


0.29 


0.04 


0.03 


8505 


18.6± 1.6 


5.77E-14 


27.07± 0.14 


-6.36 


47UMa 


GO.O V 


14.0 


5.10 


0.56 


74.0 


0.00 


2.11 


2.60 


6196 


2.6± 0.7 


1.07E-14 


26.40± 0.21 


-7.34 


51 Peg 


G5.0 V 


147 


5.49 


0.67 


37.0 


0.20 


0.05 


0.47 


25299 


0.4± 0.2 


1.70E-15 


26.65± 0.18 


-7.00 


T Boo 


M8.0 V 


15.0 


4.50 


0.48 


3.3 


0.28 


0.05 


3.90 


38251 


1252.1± 57 


3.21E-12 


28.94± 0.12 


-6.14 


HD 160691 


G3.0 IV-V 


15.3 


5.15 


0.70 


- 


0.28 


0.09 


0.04 


7046 


3.6± 1.2 


1.06E-14 


26.47± 0.16 


-7.36 


HD 190360 


G6.0 IV 


15.9 


571 


0.73 


- 


0.24 


0.13 


0.06 


2888 


2.2± 1.4 


5.23E-15 


26.20± 0.21 


-7.45 


HD 99492 


F7.0V 


18.0 


1.51 


1.01 


45.0 


0.36 


0.12 


0.11 


19928 


7.1±0.6 


2.44E-14 


26.98± 0.15 


-5.97 


14 Her 


KO.O V 


18.1 


6.67 


0.90 


- 


0.43 


2.77 


4.64 


5532 


14.6± 2.9 


3.25E-14 


27.11± 0.14 


-6.33 


HD 154345 


G8.0 V 


18.1 


6.74 


0.76 


- 


-0.11 


4.19 


0.95 


3845 


18.6± 2.4 


5.46E-14 


27.33± 0.16 


-6.03 


HD 27442 


K2.0 III 


18.1 


4.44 


1.08 




0.20 


1.18 


1.28 


4636 


37± 1.3 


1.23E-14 


26.68± 0.18 


-7.72 


ySPic 


A6.0 V 


19.3 


3.86 


0.17 


0.7 


0.00 


8.00 


8.00 


54896 


0.2± o.r 


6.00E-16 


25.40± 0.15 


-9.09 


HD 189733 


K1.5 


19.3 


1.61 


0.93 


13.4 


-0.03 


0.03 


1.13 


36271 


110.3± 1.8 


4.11E-13 


28.26± 0.12 


-4.84 


HD 217107 


G8.0 IV 


197 


6.18 


0.72 


37.0 


0.37 


0.07 


1.33 


5576 


<6.0 


< 1.55E-14 


< 26.86 


< -6.79 


HD 195019 


G3.0 IV-V 


20.0 


6.91 


0.64 


22.0 


0.08 


0.14 


3.70 


8333 


2.8± 0.8 


6.44E-15 


26.49± 0.17 


-6.86 


16 Cyg B 


G2.5 V 


21.4 


6.20 


0.66 


31.0 


0.08 


1.68 


1.68 


10768 


< 1.6 


< 5.42E-15 


< 26.47 


< -7.22 


HD 164922 


KO.O V 


21.9 


7.01 


0.80 




0.17 


2.11 


0.36 


6955 


< 3.5 


< 1.21E-14 


< 26.84 


< -6.59 


HD 4308 


GO.O V 


21.9 


6.54 


0.65 




-0.31 


0.11 


0.05 


7837 


2.1±07 


7.89E-15 


26.66± 0.19 


-6.02 


HD 114783 


KO.O 


22.0 


7.57 


0.93 




0.33 


1.20 


0.99 


3583 


2.1± 1.5 


6.72E-15 


26.59± 0.19 


-6.66 


HD 216437 


G4.0 IV-V 


26.5 


6.06 


0.63 




0.00 


2.70 


2.10 


3329 


8.2± 1.9 


1.89E-14 


27.20± 0.18 


-6.73 


HD 20367 


GO.O 


27.0 


6.41 


0.52 




0.10 


1.25 


1.07 


8861 


1404.8± 12.6 


2.76E-12 


29.38± 0.12 


-4.40 


HD 114386 


K3.0 V 


28.0 


8.80 


0.90 




0.00 


1.62 


0.99 


3601 


27± 1.2 


7.19E-15 


26.83±0.21 


-6.13 


HD 52265 


KO.O III 


28.0 


6.30 


0.54 




0.11 


0.49 


1.13 


6954 


5.6± 1.0 


1.82E-14 


27.23±0.17 


-6.92 


HD 75289 


K3.0 V 


28.9 


6.35 


0.58 


16.0 


0.29 


0.05 


0.42 


6681 


3.0± 0.7 


1.21E-14 


27.09± 0.20 


-6.75 


HD 93083 


K2.0 V 


28.9 


8.33 


0.94 


48.0 


0.15 


0.48 


0.37 


7789 


7.4± 1.3 


1.67E-14 


27.22± 0.16 


-6.79 


HD 102195 


KO.O V 


29.0 


8.06 


0.83 


12.0 


0.05 


0.05 


0.45 


13043 


145.9± 3.4 


2.87E-13 


28.46± 0.12 


-4.81 


HD 111232 


G8.0 V 


29.0 


7.61 


0.68 


30.7 


-0.36 


1.97 


6.80 


6996 


< 3.2 


<9.72E-15 


< 26.99 


< -6.41 


HD 70642 


GO.O V 


29.0 


7.18 


0.71 




0.16 


3.30 


2.00 


10935 


3.0± 0.7 


6.68E-15 


26.83±0.17 


-8.08 


HD 130322 


KO.O V 


30.0 


8.05 


0.78 


8.7 


-0.02 


0.09 


1.08 


4194 


167± 2.2 


3.87E-14 


27.62± 0.16 


-5.66 



Table 4. Stellar and planetary parameters of planet-bearing stars within 30 pc, as observed by XMM-Newton. Stellar and pl anetary 
para meters taken from www.exoplanet.eu, bolometric luminosities calculated from niy with bolometric corrections from iFloweii 
0996). 

° PN, 0.2-2 keV 

* MOS 1 countrat e given, since P N detector suffered from pile-up for this observation 
Data taken from lHempel et al] (i2005.) . combined MOS countrate given, since PN detector was optically contaminated 



planetary magnetic field itself might affect stellar wind forma- 
tion and coronal densities, as a recent study jCohen et al.ll2009l) 
suggests. 

Indic ations for SPI signatur es in stellar chromospheres were 
found bv lShkoInik et alJ (l2005h for two out of 13 stars, namely 
HD 179949 and v And, both stars with Hot Jupiters. The spectra 
of those stars showed periodic peaks in the Ca 11 H and K line 
emissions, common chromospheric activity indicators. The am- 
plitude of the variation was 2.5% for HD 179949 and 0.7% for 
V And in the K line flux compared to a mean spectrum of the re- 
spective star and appeared once per orbital period of the planet in 
several years. However, the peaks changed to a once-per-stellar- 
rotation cycle in other years, suggesting an "On/OfF'-behavior 
of chromospheric SPI. The fact that those peaks appeared only 



once per orbital period points towards magnetic SPI, since in 
a tidal SPI scenario one would expect two peaks in that time, 
which is not backed up by the Shkolnik data. 

A first statistical study on possible X- ray flux enhancement s 
due to Hot Jupiters has been conducted bv lKashvap et alJ(|2008l) . 
They claim to have found strong evidence that stars with Hot 
Jupiters are on average more X-ray active than stars with dis- 
tant planets. Their study uses main-sequence planet-bearing stars 
which were known at the time of writing, but the X-ray detection 
rate among these stars was only approximately one third, so that 
the authors had to include a large number of upper limits in their 
analysis and used Monte Carlo simulations on the X-ray lumi- 
nosities of their sample. Their analysis suggests that stars with 
planets closer than 0.15 AU have on average four times higher 
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Stai- 


Type 


Dist. 


my 


B-V 


P. 


[Fe/H] 




M„i 


logix 








(pc) 






(d) 




(AU) 


(My) 


0.1-2.4 keV 


GJ832 


Ml. 5 


4.9 


8.67 


1.46 




-0.31 


3.40 


0.64 


26.77± 0.21 


-9.09 


GL581 


M3.0 


6.3 


10.55 


1.60 


84.0 


-0.33 


0.04 


0.05 


< 26.89 


< -4.57 


Fomalhaut 


A3.0 V 


7.7 


1.16 


0.09 




0.00 


115.00 


3.00 


< 25.90 


< -8.88 


GJ 849 


M3.5 


8.8 


10.42 


1.52 




0.00 


2.35 


0.82 


27.25± 0.26 


-6.73 


HD 285968 


M2.5 V 


9.4 


9.97 


1.51 


38.9 


-0.10 


0.07 


0.03 


27.48± 0.28 


-6.13 


GJ436 


M2.5 


10.2 


10.68 


1.52 


45.0 


-0.32 


0.03 


0.07 


27.16± 0.34 


-6.36 


HD 3651 


KO.O V 


11.0 


5.80 


0.92 




0.05 


0.28 


0.20 


27.25± 0.23 


-5.66 


HD 69830 


KO.O V 


12.6 


5.95 


0.79 




-0.05 


0.08 


0.03 


27.47± 0.30 


-4.84 


HD 40307 


K2.5 V 


12.8 


7.17 


0.93 




-0.31 


0.05 


0.01 


26.99± 0.28 


-7.36 


HD 147513 


G3.0 V 


12.9 


5.37 


0.60 




-0.03 


1.26 


1.00 


29.01 ± 0.16 


-3.80 


V And 


F8.0 V 


13.5 


4.09 


0.54 


12.0 


0.09 


0.06 


0.69 


28.11± 0.22 


-6.86 


y Cep 


K2.0 V 


13.8 


3.22 


1.03 




0.00 


2.04 


1.60 


26.96± 0.20 


-7.34 


HR810 


GO.O V 


15.5 


5.40 


0.57 


7.0 


0.25 


0.91 


1.94 


28.74± 0.21 


-7.45 


HD 128311 


KO.O 


16.6 


7.51 


0.99 


11.5 


0.08 


1.10 


2.18 


28.52± 0.21 


-4.57 


HD 7924 


KO.O V 


16.8 


7.19 


0.82 




-0.15 


0.06 


0.03 


27.45± 0.29 


-6.75 


HD 10647 


F8.0 V 


17.3 


5.52 


0.53 




-0.03 


2.10 


0.91 


28.21± 0.17 


-4.88 


pCrB 


GO.O V 


17.4 


5.40 


0.61 


19.0 


-0.24 


0.22 


1.04 


< 27.69 


<-6.13 


GJ3021 


G6.0 V 


17.6 


6.59 


0.75 




0.20 


0.49 


3.32 


29.02± 0.21 


-7.00 


HD 87833 


KO.O V 


18.1 


7.56 


0.97 




0.09 


3.60 


1.78 


27.58± 0.20 


-6.79 


HD 192263 


K2.0 V 


19.9 


7.79 


0.94 


27.0 


-0.20 


0.15 


0.72 


28.03± 0.35 


-4.81 


HD 39091 


Gl.OIV 


20.5 


5.67 


0.58 




0.09 


3.29 


10.35 


27.33± 0.20 


-6.03 


HD 142 


Gl.OIV 


20.6 


5.70 


0.52 




0.04 


0.98 


1.00 


< 28.20 


< -5.63 


HD 33564 


F6.0 V 


21.0 


5.08 


0.45 




-0.12 


1.10 


9.10 


27.84± 0.30 


-6.66 


HD 210277 


GO.O V 


21.3 


6.63 


0.71 




0.19 


1.10 


1.23 


< 27.85 


< -5.68 


70 Vii 


G4.0 V 


22.0 


5.00 


0.69 


31.0 


-0.03 


0.48 


7.44 


27.42± 0.28 


-6.33 


HD 19994 


F8.0 V 


22.4 


5.07 


0.57 




0.23 


1.30 


2.00 


28.16± 0.28 


-6.41 


HD 134987 


G5.0 V 


25.0 


6.45 


0.70 




0.23 


0.78 


1.58 


< 27.99 


< -5.75 


HD 16417 


Gl.O V 


25.5 


5.78 


0.67 


_ 


0.19 


0.14 


0.07 


< 28.28 


< -5.73 


HD 60532 


F6.0 IV-V 


25.7 


4.45 


0.48 




-0.42 


0.76 


3.15 


< 26.98 


< -7.53 


HD 181433 


K3.0 IV 


26.1 


8.38 


1.04 




0.33 


0.08 


0.02 


< 27.08 


< -6.05 


HD 30562 


F8.0 V 


26.5 


5.77 


0.63 




0.24 


2.30 


1.29 


< 26.97 


< -7.07 


HD 179949 


F8.0 V 


27.0 


6.25 


0.50 


9.0 


0.22 


0.05 


0.95 


28.61± 0.25 


-4.40 


HD 150706 


GO.O 


27.2 


7.03 


0.57 




-0.13 


0.82 


1.00 


28.88± 0.19 


-5.04 


HD 82943 


GO.O V 


27.5 


6.54 


0.62 




0.27 


0.75 


2.01 


< 28.01 


< -5.75 



Table 5. Stellar and planetary parameters of planet-bearing stars within 30 pc, as observed by ROSAT. Stellar and p lanetary nararn - 
eters taken from www.exopla net.eu, bolometric lum inosities calculated from my with bolometric corrections from lFlowerl(ll996h . 
X-ray luminosities taken from'K ashvap et alj (|2008|) . except for upper Umits for HD 1641 7, HD 30562, HD 181433 and HD 60532, 
which were calculated from original data. 



X-ray luminosity than stars with planets at distances larger than 
1 .5 AU. They try to account for selection effects by regarding the 
trend of Lx/Lboi with a pi as selection-induced and the remaining 
trend in Lx with Opi as planet-induced, which leads to remaining 
Lx-ratio of stars with close-in and far-out planets of x 2, with an 
overlap of the (simulated) Ljf-distributions at Icr level. 

We do not see a significant difference of Lx-distributions in 
dependence on a pi in our sample as shown in Fig.|4] There is also 
no significant trend of LxjLhoi with a pi evident in our data. We 
do not try to correct artificially for selection effects, since these 
are various and interdependent: since stellar activity masks the 
planet-induced radial-velocity signal, small far-out planets are 
more easily detected around very inactive stars. Similarly, those 
planets are also easier to detect around low-mass and therefore 
late-type stars, but very late-type stars have again higher LxlLi,„i 
values than earlier-type and therefore heavier stars. A quanti- 
tative estimate of activity-related selection effects is therefore 
extremely difficult. But as shown in section l6n there is no sig- 
nificant additional effect on activity visible in our data which can 
be attributed to the influence of massive close-in planets. 

There has also been an effort to measu re coronal SPI fo r 
an individual target; the star HD 179949 bv ISaar et aP (l2008h . 
which showed the largest SPI signatures in chromospheric data 
so far. The star's X-ray flux was measured several times during 
May 2005; in September 2005, the star was in an "On"-state of 
SPI as seen in the chromosphere. The measured X-ray fluxes 
above 0.3 keV vary by +15%, a typical level also for intrinsic 



stellar X-ray variability. When folded with the orbital period, 
the profile of X-ray variability does not match the one seen in 
chromospheric data very well; interpretations of variability with 
the stellar rotation period are also possible. To see how chromo- 
spheric and coronal variability compare with each other for this 
star, we can do a rough estimate: a variation of 2.5% in the Can 
K line compared to the mean stellar spectrum should translate 
into something of the same relative order of magnitude for the 
Mount Wilson S index. The S index variation should even be a 
bit smaller than 2.5%, since the S index averages over the H and 
K line, and the H line is less sensitive to activity effects. If we 
compare this to stars with known activity cycles such as 61 Cyg, 
one finds there (iHempebnann et al ] 12006) that the S index dur- 
ing one stellar activity cycle varies by +15%, while the X-ray 
flux in the 0.2 - 2.0 keV band varies by +40%. A similar ratio 
of X-ray to Ca ii fluxes yields for HD 1 79949 an expected X-ray 
variation of ^ 7%, less than a typical intrinsic variability level 
for late-type main-sequence stars. However, this ratio between 
X-ray and Ca ii fluxes should only apply if activity enhancement 
via SPI works via similar mechanisms as normal stellar activity 
does, which is not necessarily the case given the possibility of 
Jupiter-Io-like interaction scenarios. 

A basic quant itative scenario for coronal magnetic SPI has 
been proposed by ICohen et al.l (|2009'), suggesting that the pres- 
ence of a close-in planetary magnetosphere hinders the expan- 
sion of the stellar magnetic field and the acceleration of the stel- 
lar wind, causing a higher plasma density in a coronal "hot spot". 
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In their model, the hot spot leads to variations of the X-ray flux 
of 30% when rotating in and out of view, and compared to a 
setting without a planet, the overall X-ray flux is enhanced by a 
factor of a: 1.5 for a stellar dipole field and a factor of w 15 for 
a stellar magnetic field like the Sun's in an activity maximum. If 
we compare this to our sample, we can conclude that an activ- 
ity enhancement of more than one order of magnitude is not a 
common effect in stars with Hot Jupiters, since such an activity 
overshoot for stars with Hot Jupiters would yield significantly 
different results in our sample's trend of Lx with log(fl^^^' x Mpi). 

In summary, we can conclude from our analysis that there 
is no major average activity enhancement in the corona of stars 
which is induced by their planets. Any trends seen in our sample 
seem to be dominated by selection effects. 

6.3. Promising individual targets 

SPI has been claimed for a handful of targets in chromospheric 
emissions at selected times. Even if there has been no stringent 
detection of corresponding coronal SPI yet for these stars, si- 
multaneous observations in the optical and X-ray regime yield 
precious insight into the interplay of stellar and planetary mag- 
netic fields. The most promising candidates for such coordi- 
nated searches are stars with close-in and massive planets which 
are rather X-ray bright to allow feasible observation plans. This 
identifies the stars HD 102195, HD 130322, HD 189733, t Boo, 
V And, HD 179949 as promising targets. Four of these, namely 
HD 189733, t Boo, v And and HD 179949, have been in vesti- 
gated for chromospheric SPI before CShkolnik et al.ll2005h . Out 
of these, u And and HD 179949 did show chromospheric ac- 
tivity enhancement in phase with the planetary orbit; the other 
two stars showed indications for increased variability of chro- 
mosp heric activity with Port dWalker et al.ll2008l:IShkolnik et alj 
l2005h . The remaining two stars, HD 102195 and HD 130322, 
have not been analyzed in detail for chromospheric SPI yet, but 
might as well be interesting targets for combined chromospheric 
and coronal SPI searches. 



5. Coronal SPI might still be observable for some individual 
promising targets, preferably in coordinated observations of 
the targets' coronae and chromospheres. 
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7. Conclusions 

We analyzed a sample of all known planet-bearing stars in the 
solar neighborhood for X-ray activity and possible manifesta- 
tions of coronal Star-Planet-Interactions (SPI) in dependence of 
the planetary parameters mass and semimajor axis. Our main re- 
sults are summarized as follows: 

1 . In our sample of 72 stars, there are no significant correla- 
tions of the activity indicator Lx jLj,,,! with planetary mass or 
semimajor axis. 

2. However, we do find a correlation of the X-ray luminosity 
with the product of planetary mass and reciprocal semimajor 
axis. Massive close-in planets are often found around X-ray 
brighter stars. 

3. This dependence can be ascribed to selection effects: The ra- 
dial velocity method for planet detections favors small and 
far-out planets to be detected around low-activity. X-ray dim 
stars. Additionally, if SPI induced an excess in Lx without 
changing the Lx/Lhoi ratio, SPI would need to cause ex- 
tremely high energy input in Lhoi, leading to unrealistically 
short decay times for the planetary orbit. 

4. There is no additional effect detectable in Lx which could be 
attributed to coronal manifestations of SPI. 



